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1 . SUMMARY 

Compressor stall phenomena are analyzed from a nonlinear control 
theory viewpoint, based on bifurcation-catastrophe techniques. This new 
approach appears promising and offers insight into such well known 
compressor instability problems as surge and rotating stall; furthermore it 
suggests strategies for recovery from stall. More specifically, three 
Interlocking dynamic nonlinear state space models are developed. The 

models, which are based on Wenzel, Bruton (1982), represent compressor 
operation under normal unstalled conditions (US), and stalled compressor 
operation under normal and reverse flow conditions (PS and RS, 
respectively). It is shown that the problem of rotating stall can be 

viewed as an induced bifurcation of the solution of the (US) model. Both 
equilibrium and bifurcation loci are generated in the state-parameter 
space, based on parameter continuation techniques combined with 

catastrophe/singularity theory. Hysteresis effects are shown to exist in 
the stall/recovery process. Surge cycles (of about 50 cycles/sec) are 
observed for some critical parameter values. It is shown that the oscilla- 
tory behavior is due to development of limit cycles, generated by Hopf 
bifurcation of solutions. More specifically it is observed that at certain 
critical values of parameters a family of stable limit cycles with growing 
and then diminishing amplitudes is generated, then giving rise to an 
unstable family of limit cycles. This unstable family in turn bifurcates 
into other unstable families (as indicated in the diagrams). The quan- 
titative approach adopted for generation and continuation of limit cycles 
employs a number of notions from topology, differential calculus, and 
singularity theory. (Unstable limit cycles can not be detected by simula- 
tion techniques, though their presence vastly influences the system 
dynamics.) To further illustrate the usefulness of the methodology some 
partial computation of domains of attraction of (US) and (PS) equilibria is 
carried out, though this objective is clearly beyond the scope of this 
investigation. This provides a global picture, suitable to address the 
non-steady state 'off-design' instability conditions which are present 
before the onset of stall. 

The domains of attraction of equilibria help define the 'position' of 
system trajectories in a global sense, and in this way provide a natural 
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classification of trajectories with respect to their stall/recovery charac- 
teristics. Three cypes of trajectories are defined according to the 'se- 
verity' of their stall/recovery action: 1st stage recoverable, 2nd stage 
recoverable, and 3rd stage recoverable. Finally, it is argued that the use 
of corrected flow (e.g. a scalar output variable), as the sole indicator of 
stall behavior may lead to incorrect conclusions with respect to recovera- 
bility of trajectories. In particular it is argued that neither the 
corrected flow rate nor the turbine area (or any other single scalar 
variable) can provide an unambiguous measure of stability. The use of the 
domain of attraction of some normal operating point(s) is proposed as a 
useful mechanism for stall detection, and recovery. By varying the parame- 
ters over their feasible range, the size and shape of the domains of 
attraction change, leading to the determination of those trajectories which 
may be recovered by appropriate control actions. Furthermore this approach 
provides a natural framework for carrying out parameter sensitivity analy- 
sis yielding the impact of various parameter changes on recovery from 
stall. Some sample sensitivity computations are also offered. 
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2. INTRODUCTION 

Stabilicy of aircraft engines over their operating envelopes is a pri- 
mary design consideration. The degree of stability varies due to a number 
of factors such as engine operating condition, inlet flow conditions, age, 
and control tolerance and is generally a trade-off with performance; that 
is, for a given engine, increases in performance are attained at the 
expense of reduced stability margin. As a result, gas turbine engines have 
experienced stalls or surges throughout their history as aircraft 
power plants. 

The stall problem has been receiving increasing attention in recent 
years. It was the subject of a recent NASA Lewis Workshop (1983). The 
major aircraft engine manufacturers such as GE and Pratt & Whitney are 
developing complex simulation models, whereas research at universities 
(MIT, Cornell) is aimed at developing a more basic understanding of the 
problem. 

Most stalls result in momentary disturbances to engine operation and 
recover without any active intervention. Some stalls are of more serious 
nature and require engine shutdown to restore normal engine operation. 
Such stalls are typically characterized by low engine speed and high tur- 
bine temperatures together with low thrust; these stalls are referred to as 
nonrecoverable , stagnation, or hung stalls. Prolonged operation in nonre- 
coverable state can result in turbine damage due to over-temperature; in 
this state the engine is not responsive to scheduled inputs. Indeed, as 
expressed in Stetson (1983), "one of the most critical functional problems 
that a high technology turbine engine encounters is nonrecoverable stall." 

In recent years, the incidence rate of nonrecoverable stall has 
increased by a great deal. This is mainly due to the increased use of 
augmented, mixed-flow, turbofan engines in modern aircraft due to their 
wide performance range, compact size, and flexibility. In these aircraft 
the engine is a highly interactive system. This interaction between com- 
ponents can directly affect the stall/surge characteristics through compli- 
cated feedback paths (Patterson (1983)). 
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In Che foLiowing paragraphs we discuss che above issues in some more 
detail. 

Typically, multimission aircraft experience a wide range of Mach num- 
bers, Reynolds- numbers, angles of attack., side slip angles and dynamic 
motions in longitudinal and lateral planes. However, aircraft, including 
the propulsion system, have historically been designed for steady con- 
ditions with minimum attention to the effect of distorted flow. Actual 
operation of the aircraft is in an environment which is unsteady, and in 
which the propulsion system experiences distorted flows. 

The propulsion system, consisting of the gas turbine engine, is par- 
ticularly sensitive to these "off design" operating conditions. In parti- 
cular, it is the compressor which has a narrow operating range. If one 
exceeds the limits of the operating range, the compressor operation becomes 
unstable. Such instabilities may affect the flow in the entire engine. 
The combustor or afterburner may discontinue burning, resulting in loss of 
thrust for the aircraft. 

The result of many of these instabilities is only a temporary disrup- 
tion (less than a second, say) of the unstalled operation of the engine. 
In some cases, however, the consequence of the initial instability is a 
situation in which the engine will not recover to unstalled operation and 
will have to be shut down and restarted in flight, and this has an obvious 
implication for the effectiveness of the aircraft. 

If such a non-recoverable stall occurs the sequence of events is that 
the turbine inlet temperature will rise and the engine speed will fall. 
The rate at which the former occurs is determined from the interactions of 
the fuel control and the compression, burner and turbine system. The drop 
in speed arises due to the mismatch of turbine and compressor work, because 
of the very low compressor efficiencies in this regime. The danger from 
non-recoverable stall is that the combustor exit temperature (turbine inlet 
temperature) will exceed the allowable limits for the turbine and/or that 
the rotor speed will fall below the level that the gas turbine power cycle 
is self-sustaining. If so, it is generally necessary to shut the engine 
down and cool it before restarting. 
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As stated this problem has received a great deal of attention in 
recent years, and there has been much work on the topic, at levels ranging 
from examination of the detailed fluid mechanics of unsteady boundary 
layers on compressor blades to analysis and experiments which attempted to 
consider a description of the overall system. These latter results have 
proved to be of significant help in understanding the overall phenomenon. 

A basic analysis of this type was presented in Greitzer (1976) in 
which the nonlinear behavior of a model compression system was examined. 
Experiments were also carried out to validate system modeling. Since that 
study, system models have also been carried out by many other investiga- 
tors. 

At NASA Lewis, Wenzel and Bruton (1982) have developed a lumped para- 
meter simulation model of the TF34 engine to study the nonrecoverable stall 
phenomenon. This study, in fact, forms the basis of the models presented 
in this report. Other relevant NASA reports include: Szuch and Seldner 
(1975), Szuch, Seldner and Cwynar (1977), Szuch (1978); and Seldner, 
Mihaloew and Blaha (1972). A major theoretical contributor is Frank Moore 
(1983). Another contribution to the subject is the classic work of Emmons 
et al. (1955), which appears to have laid the foundation for much of the 
current research in the compressor stall area. 

Though, as the bibliography indicates, there has been considerable 
interest in the subject of stall, the phenomenon is still not well 
understood, especially with regard to the system behavior subsequent to the 
onset of instability, i.e. , the post stall response of the system. This is 
true not only as far as tests and/or experiments are concerned, but also as 
far as modeling and analytical understanding are concerned. This is in 
part due to the fact that .the post-stall transients involve very compli- 
cated nonlinear interactions and feedback mechanisms — a heavily mathema- 
tical subject. In particular, up until the present there has not been an 
attempt to bring to bear on this problem some of the modern mathematical 
techniques employed in nonlinear control theory. This report presents such 
an analysis. 
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This work addresses the application of bifurcation-catastrophe theory 
to the problem of nonlinear compression system stability. It is thus an 
attack on the problem from a different viewpoint than existing work in this 
area. Hence it provides a complementary perspective for analyzing the 
strongly nonlinear dynamics of compression systems. 

The report is organized as follows. Chapter 3 presents some 
background information on bifurcation theory used in subsequent chapters. 
Chapter 4 presents three state space models of compressor dynamics: one 
for the unstalled compressor, and two for the stalled compressor — 
corresponding to normal flow and reverse flow conditions. (Appendices B, 
C, and D present a derivation of the models based on the modified 
Wenzel-Bruton (1982) results.) In Chapter 5, generation of equilibrium 
surfaces is discussed. Chapter 6 is devoted to the analysis of 
bifurcation-catastrophe points of the models; Hopf bifurcation points are 
discussed. Issues regarding stall recovery are discussed in Chapter 7; 
three types of trajectories corresponding to their stall/recovery behavior 
are defined. Finally, Chapter 8 contains some concluding remarks. 
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3. PERTINENT BACKGROUND IN BIFURCATION TECHNIQUES l 


3.1 Introduction 


While linear control theory is well developed by now, the same asser- 
tion cannot be made about the nonlinear theory. The problem is even more 
frustrating if one is interested in obtaining computable results in par- 
ticular. This is due to the fact that nonlinear theory tends to be more 
qualitative . 

One of the most exciting developments in the field of nonlinear stabi- 
lity theory has been that of bifurcation theory — pioneered by the work of 
Poincare, and later refined by Andronov (1937, 1966), Hopf (1942), and 
others. The Importance of this theory is in its practical applicability. 
For example, in elastic stability, bifurcation is called buckling, and in 
hydrodynamic stability bifurcating solutions are referred to as tran- 
sitions. It has numerous applications in many other physical systems. 

Another important feature of the theory is in its ability to provide a 
global picture. That is, it gives information about all of the "turning 
points" of the solutions. This is in contrast to linearization procedures, 
in which only local stability is characterized. In fact, in a certain 
sense, bifurcation theory takes over when linearization fails. 

3 .2 Hyperbolic Fixed Points and Linearization 

We shall use the above remarks to motivate our discussion. Consider 
the nonlinear differential equation 


dx(t) 

dt 


f(x,p) 


(3.1) 


*This chapter is based on SSI's Progress Report (Jan. I - May 1, 1984) to 
NASA Lewis. 
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in which x X C R n 

u E U C R m . 

We shall refer Co x as the state vector and u as the control vector (or 
sometimes as the parameter). For the purposes of illustration let y = 0, 
i.e., consider the system 

= f(x), x e X C R n (3.1)' 


Assume that the point xq£ R n is a stationary point of f in (3.1)': 


0 = f ( xq ) 


(3.2) 


and let the Jacobian matrix of partial derivates of f at xq be denoted as 
Vf(x 0 ). 


Consider the linear system 


(3.3) 


where 

z = x 0 + e,| 

The major question to be addressed HT~what we 
tions of (3.1)' based on our knowledge of (3.3)? 
interest since the solutions of (3.3) are given by: 



(3.4) 

can say about the solu- 
This question is of 


z( t) 


g tVf(x 0 ) 


z(0) 


(3.5) 


To relate (3.5) and solutions of (3.1)' let <(> t (x) denote the flow generated 
by f, i.e., a smooth function defined for all x e X and t e (a,b) C R such 
that 

= f(<H( x )) (3.6) 
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for all x : X and c z (a,b). We now state a fundamental result in nonli- 
near system theory which provides the answer. 

Result 1 (Hartman ( 1964 )) 

Provided that 7f (xg) has no purely imaginary or zero eigenvalues, there 
is a continuous mapping with a continuous inverse (i.e., a homeoraorphism) , 
say, g, defined on some neighborhood of xg in R n locally taking trajec- 
tories of the nonlinear flow <j> t of (3.1)' to those of the linear flow 
e f( x o) c 4 The homeomorphism preserves the direction of the trajectories. 

We will not provide the proof of this result here. The interpretation 
is that roughly speaking, as long as the linearized system has no purely 
imaginary or zero eigenvalues, it provides a good local picture of the 
nonlinear system. (It is possible to strengthen this result in the sense 
that the linear and nonlinear system solutions may be related via a dif- 
feomorphism. But this requires certain "resonance" conditions among the 
eigenvalues. We shall not discuss the matter here.) The important point 
for our purposes is that when Vf(xg) has no eigenvalues with zero real 
parts it is acceptable to linearize. A stationary point xg at which Vf(xg) 
has the above property is known as being hyperbolic. (But of course bifur- 
cating solutions occur precisely when this condition is not met!) The 
following result will further characterize hyperbolic fixed points (or sta- 
tionary points). But first we need the following definitions. Define the 
local stable and unstable manifolds of xq , 

M s (xg) = {x e x such that 4>t( x ) + xg is t + ® 
and t ( x ) e X for all t > 0} 

M u (xg) = (x e X such that <j> t (x) xg as t ♦ -« 
and $ t (x) e X for all t < 0} 

Let E s . and E u denote the eigenspaces associated with stable and 
unstable eigenvalues. They are respectively said to be stable and unstable 
eigenspaces . 
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Result 2 (Hyperbolic Fixed Points) 

Let xq be a hyperbolic fixed point of (3.1)’. Then there exist local 
stable and unstable manifolds M 3 (xq), M u (xq), which are, respectively, of 
the same dimensions d a and d u as those of the eigenspaces E s and E u of the 
linear system (3.3). Furthermore, the manifolds M 3 (xq) and M^xq) are 
tangent to the eigenspaces E 3 and E u , respectively. The manifolds are as 
smooth as the function f. 



Figure 3.1 Manifolds M 3 , and M°; Eigenspaces E 3 and E u . 
3.3 Center Manifolds 


As mentioned previously, in bifurcation problems we are specifically 
interested in equilibria at which we have purely imaginary -eigenvalues. On 
the complex plane these eigenvalues are located on the imaginary axis and 
separate the stable eigenvalues (e.g., those with negative real parts) from 
the unstable ones, (e.g., those with positive parts). Hence they are 
located in the "center". Thus in the linear analysis attention is focused 
on this center subspace associated with the purely imaginary eigenvalues. 
In nonlinear problems we are interested in center manifolds. The other two 
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manifolds are referred to as stable and unstable manifolds. 8y focusing on 
this manifold we will in fact reduce the dimension of the problem in a 
systematic manner. 

The center manifold concept isolates the complicated asymptotic beha- 
vior by pinning down an invariant manifold tangent to the subspace spanned 
by the generalized eigenspace of eigenvalues on the imaginary axis. But 
one problem here is that center manifolds are in general not unique. This 
is in contrast to the stable manifolds. 

The above discussion is made more precise in the following statement. 

Result 3 (Center Manifolds) 

Let f be k times differentiable (denoted by C k ) with f(0) = 0. (This 
may be accomplished by a shift of origin.) Thus the Jacobian in (3.3) 
becomes the square matrix Vf(0). Split the spectrum of Vf(o) into o s , o c , 
o u such that 

Re X > 0 if X £ o u 
Re X = 0 if X s o c 
Re X < 0 if X e a s 

with E s , E u , and E c denoting the respective eigenspaces of o s , o u , and o c . 
Then there exist stable and unstable invariant manifolds M s and M 11 
tangent to E s and E u at zero. Also there exists a C^-l center manifold M c 
tangent to E c at zero. Both the stable and the unstable manifolds are uni- 
que. But in general M c may not be unique. 

The following example, due to Kelley (1967) illustrates this point. 
Consider the differential equation 



x 2 = -x 2 


It is easily verified that the solutions to this system are 
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Eliminating t yields 


x L (t) 

X2(t) 


ci 





-_L_ J_ 

. v C 1 X 1 

x 2v x l) = C 2 e e 1 



Figure 3.2 Center Manifold Is Not Unique 


It is clear that for x > 0 the only solution approaching the origin is 
x = 0. But for x < 0, any solution with zero slope at the origin is accep- 
table. For example, the entire x^ axis is one such solution, and in fact 
the only analytic one. This is illustrated in Figure 3.2. (Note that a 
center manifold is an invariant manifold tangent to the center eigenspace. 
Hence this tangency condition at zero can be satisfied by many center mani- 
folds: take any solution curve with xj<0 pieced together with the positive 
half of the x^ axis; the resulting curve will be a center manifold.) 
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Figure 3.3 Stable, Unstable, and Center Manifolds 


3 .4 Hopf Bifurcations 

The simple equilibria that we have discussed so far are examples of 
point equilibria. It is possible, for example, to have situations in which 
the state of a system hovers around in the state space without ever 
reaching any stationary point. One such situation arises in the following 
example . 


■— ■ 


' ' 


— r 

*1 

3 

a l “ a 2 


X 1 

CM 

•X 

a 

02 + a ^ 


x 2 

U - 


__ 


__ 


The solutions to (3.8) are of the form 


(3.8) 


14 



x^(t) 


X2(t) 


= e a l C 


cos a2 t -sin <*2 t 


sin a2 t cos 02 t 


= Scientific Systems 


*l(0) 


* 2 (0) 


(3.9) 


For > q t ^ e trajectories spiral away from the origin; for < 0 the 
opposite occurs. Bur for = 0 the solutions become periodic — never 
expanding or contracting. The Hopf theory is a body of results which 

generalize this observation. Note that the eigenvalues of the matrix 


a l -o 2 
o 2 O! 


are ± ja2 and become purely imaginary for aj^ = 0. 
Let's slightly generalize the form of (3.8): 


*1 * (dai + a(x^ + x|)) x^ - (02 + caj + b(x^ + x|))x2 

*2 = ( a 2 + ca l + b( x | + x |)*l + ( dx l + a ^ x i + x 2^ x 2 


or expressing (3.10) in polar coordinates 


r = (doj + ar^)r 


0 = (02 + caj + b r^) 


(3.10) 


(3.11) 


where 
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r = /x2 + 

* 

0 = tan l — 
x 2 


In equation (3.11) the radius r is independent of 0 and therefore r=0 
in (3.11) yields the steady state trajectories — which are circles of 
constant r. 

In a fundamental sense, equation (3.11) captures the essence of the 
Hopf bifurcation in that the addition of higher order terms will not make 
a qualitative difference in the behavior of solutions near the origin: 

Result 4 (Hopf Bifurcation) 

Consider the nonlinear control system (3.1), and assume that it has a 
steady state solution xq corresponding to p =■ mq with the property (PI): 
7 x f( x O» y o) 1133 a simple pair of purely imaginary eigenvalues and no other 
eigenvalues with zero real parts. 

Then there is a smooth curve of equilibria (x(p),pq) such that 
x(Uq) = xq . The imaginary eigenvalue pairs X(p) and X(u) vary smoothly 
with u . Furthermore, if 


•^j<Re X(u)) - d j* 0 (3.12) 

Then there is a three dimensional center manifold which passes through 
(xq.^o) in the product space R n xR. Furthermore, when (3.12) is satisfied, 
there is a -smooth coordinate system in which the Taylor series expansion of 
degree 3 on the center manifold is given by (3.10) (with the new notation 
u replacing aj). In addition, if a.* 0, there is a surface of periodic 
solutions in the center manifold tangent to the eigenspace of A(pg) (and 
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X(Uq)) which agrees to the second order with the paraboloid y = 



For a < 0, these periodic solutions are stable limit cycles — while for 
a > 0 the limit cycles are diverging. The proof of this well known result 
can be found in any standard text on nonlinear dynamic systems (c.f. 
Marsden and McCracken (1976)). 


We will see in Section 6.3 that, for some critical parameter values, 
the unstalled compressor exhibits surge cycles corresponding to Hopf bifur- 
cation points; also in that section a more quantitative discussion of Hopf 
solution points will be carried out. Mucl r of the stationary operating 
points of the compressor models will turn out to be of the Hyperbolic 
Stable or Unstable type (e.g. Section 5.5) 
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4. THREE INTERLOCKING MODELS 


4.1 Introduction 


In this chapter three separate state space models of a compressor rig 
are developed. One corresponds to the compressor operation under normal 
unstalled conditions (US); the other two are developed for the stalled 
compressor - corresponding to the normal (positive) flow and reverse flow 
conditions ((PS) and (RS), respectively). The reverse flow and normal flow 
stall conditions are modeled separately to account for the different (e.g., 
opposite) heat and mass flow conditions in the compressor. Wenzel and 
Bruton (1982) describe these conditions and provide explanations for the 
approach; their work forms the basis of the models presented here. 


The (simplified) Wenzel-Bruton model reported in Appendix A is used 
exclusively in our development of the state space models. The reader is 
assumed to be famiLiar with the related fluid mechanic terminology and 
vocabulary as described in Wenzel and Bruton (1982) or similar treatments. 
Refering to Appendix A (and Wenzel and Bruton (1982) for definitions) we 
begin with defining the state vector x as 


w 2 

(WT) 2 
(L 2 
w 3 

(wt) 3 


(4.1) 


our objective is to derive from the relations of Appendix A, three state, 
space formulations of the form 


= fi(x(t), 0(t), u(t)) (4.2) 
i = 1, 2, 3, 

in which x, 9, and u represent the state vector, the parameter vector, and 
the control vector, respectively. In addition a set of output relations of 
the form 
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Yi = gi(x» 3 » u) (4.3) 

i =1,2, 3, 


are also obtained. 

A typical element of the control vector may be the area A determining 
the flow rate downstream of the compressor. A typical parameter would be 
the fractional corrected speed PN, or some design parameter such as in- 
stall compressor performance at positive (or negative) flow KP (or KN) . 
The elements of the state vector are defined in (4.1). The first two ele- 
ments denote the mass and (thermal) energy stored upstream of the 

compressor; w 2 represents the mass flow rate through the compressor; and 
the last two elements are the downstream mass and energy variables. 

4 .2 State Space Relations for Positive Flow Stalled Dynamics (PS) 

Keeping in mind the definition of the state x as described by Equation 
4.1, Appendix B shows how the integral relations of Appendix A can be used 
to derive a state space model of the stalled compressor rig in normal 
(positive) flow condition. They are summarized below. 

Summary of State Space Systems for Positive Flow Stalled Dynamics (PS) 

. P 1 R , 

*1 * Rj * *2 ‘ ”3 <PS1) 


*2 ' * £ ^ Tl * 2 ' X2 M 

2 

*3 3 £ 83 81 *P £ 82< p N) x 2 - -J- x 5 


( PS2) 


( PS3) 


x 4 = x 3 - A — /x 5 X 4 


(PS4) 
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0.286 


X5 = r 


*2*3 

X 1 


1 + 


1.19 + 0.925 KP LSd X J.) 2 + 0.33 PN 2 - 1 

519 xix2 


n 0 + F4( pN * n 0)^ 


x 3 


/ 519 x^X2 


YAR 


x 5 


V 3 X5 x 4 

where g's are defined in Appendixes B, C, D, 

4 .3 State Space Relations for Unstalled Dynamics (US) 

Recalling from Appendix A that 


(PS5) 


Fi( 


(i>2^ 9 2 


, PN) = 


(Ilo/ITT 

f 2 ( -nrA PN > = 


p 3 * 

pT 

t 3 * 


(4.4) 


(4.5) 


«2 T 2 

we now summarize the state space model derived in Appendix C for the 
unstalled compressor. 


Summary of State Space Systems for Unstalled Dynamics (US) 


X 1 


R o 


_R 

r o v 2 


x 2 - x 3 


(Ul) 


x 2 = Y 


R 

x 3 = r 


p l 


^2 

v 2 


r qV2 


x 2 


Ti - x 3 x 2/ x l 


* F 1 S 4 


x 3 


XlX 


PN - 


l x 2 


ill 

v 3 


( 02 ) 

(03) 


. AR / 

x 4 = x 3 " yj /x 4 x 5 



(04) 

(05) 
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where 


S4 


14.7 V 2 
R /Tf 


and corrected flow 


x 3 


x l x 2 


4 .4 State Space Relations for Reverse Flow Stalled Dynamics (RS) 

The major distinction between the (PS) case and the present one is that 
in the latter the flow through the compressor reverses — thus, roughly 
speaking, the exit and inlet conditions are switched. In this case the 
compressor flow is corrected by downstream conditions. However, the 
correction to compressor speed is not switched, to preserve continuity at 
zero flow. The detailed derivation of the relations are given in Appendix 
D. Here we summarize the relations. 

Summary of Reverse Flow Stalled Dynamics (RS) 


Xi = 


P[ R 
Rq RqV 2 


x 2 - x 3 


(RSI) 


x 2 



(RS2) 


x 2 


X 3 = L 1 0.837 V 2 '\ g 9 + ^ *10 


2 

x 3 

X 4 X 5 ) 



(RS3) 


. AR / 

x 4 = x 3 ~ /x 4X5 


(RS4) 
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X 5 = Y U 3 


x 5 AR 

' v 3 


x 5 



where 


corrected flow h r = gg 


x 3 

J X4X5 


(RS5) 


Pi v 3 



(Other g's are defined in Appendix D.) 

4 .5 Typical Output Relations 

The dynamic relations (US), (PS), and (RS) describe the compressor in 
unstalled and stalled conditions. The state vector x in each case comple- 
tely summarizes the past history of the variables and can be used to 
generate future conditions under alternative input scenarios. In general, 
the state variables are not all measurable — they represent "internal" 
variables of the system. More often some other related variables, typi- 
cally represented by a combination of the state variables, are measurable 
or physically significant. They may be labled as output variables. The 
output variables y^ may be written in the general form 

yf(t) = gj.(x(t), 9(t), u(t)) 

where i varies over the set of output relations of the models (US), (PS), 
and (RS). 

In the context of our work, typical output variables are the various 
pressures and temperatures — both in absolute and relative terras. Another 
important variable that is interpretable as an output variable is the 
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Cj 2/ 9 2 

corrected flow — ? . This variable plays a crucial role in our analy- 

6 2 

sis, as will later be described. Here we cite some typical output rela- 
tions (derived in Appendices B, C, and D). 


P 2 = - 2 *2 


( 4 . 6 ) 


*2 =T7 


2! 

x i 


( 4 . 7 ) 


?3 


v7 x 5 


( 4 . 8 ) 


T 3 = — 

J V / 


ill 

x 4 


( 4 . 9 ) 


hi : = 


u>2*^ 

TT = g4 


x 3 


X 2 X 2 


(US, PS) 


( 4 . 10 ) 


i 2 /e7 x 3 

hr: = rr = g6 


2 f ^ 4 x 5 


(RS) 


( 4 . 11 ) 


* 

P3 


pj “ Fl<M,PN) 


(US) 


( 4 . 12 ) 


* 

7T = F 2 ( hi ,PN) 


(US) 


( 4 . 13 ) 


♦ = gl 


/T 


_^3 

1 X 1 x 2 


(PS) 


( 4 . 14 ) 


- 23 - 



Scientific Systems 


P 3 (gjX 3> 2 

= 1.194 + 0.925 KP + 0.33 PN 2 (PS) 

P 2 Ti x 2 


(4.13) 


12 

t 2 


1 + 


(g, X 3) 2 0.286 

1.19 + 0.925 KP + 0.33 PN 2 - 1 

519 xjx 2 


HQ + F4( pN » n 0^ 


x 3 


/519 xjx 2 


87 x 3 

PN , 

» X4X5 


(PS) (4.16) 

(RS) (4.17) 


~ = 1.194(l + 0.2556 PN 2 + 0.7746 KN (g7X3) 


p 2 


x 4 x 5 


(RS) 


(4.18) 


T* 

Tj = 1 + 


„ (g 7 X 3) 2 0.286 

1.19( 1 + 0.2556 PN 2 + 0.7746 KN ) - 1 

V5 

gl X 3 

n 0 + F^pn.hq) — 


PN / 


X 4 X 5 


(RS) (4.19) 

Despite the seemingly complicated appearance of the above functions, it is 
clear that they are all of the basic (static) output form mentioned 
earlier. 
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4.6 Linking the Models 


Having presented the three state space systems (US), (PS), and (RS), we 
now describe how they are linked. 

The criterion used for switching the models is the magnitude and direc- 
tion of the corrected flow through the compressor. For a given corrected 
speed PN, if the corrected flow falls below a certain positive threshhold 
value t 3 (PN), the compressor is assumed to enter into stall (PS). At a 
given speed PN, the requirement for recovery is for the corrected flow to 
exceed a certain threshhold t r (PN). At the time the flow reverses through 
the compressor, the (RS) dynamics are switched on; upon reentry into posi- 
tive flow region, the (PS) dynamics are switched back on. This process i3 
described graphically in Figure 4.1. 


hr > c hi 2 t r (PN) 



In closing we note that the switching of the models is based on the value 
of the corrected flow, which is an output variable. There are some 
possible criticisms of this approach — some of which we shall mention 
later. 
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5. EQUILIBRIUM LOCI 


5 . 1 Introduction 

In this chapter we will trace the equilibrium solutions of the systems 
(US) and (PS), continuously varying a parameter of interest. We will also 
show that (RS) has _no equilibria — a result which turns out to be quite 
intuitive. Various types of fixed points such as hyperbolic stable and un- 
stable points are encountered and discussed. The numerical results and 
plots are obtained by the use of our in-house software package BISTAB; a 
discussion of its features is also offered in this chapter. 

5 .2 (RS) Model Contains No Equilibria 

In order to see that the in-stall reverse flow model (RS) can not have 
any equilibria let us assume on the contrary that it does. This means, in 
particular, that the left hand side of Equation (RS4) must equal Eero, 
1 . e . , 

0 = K 3 - ^* / X 4 X 5 

with X3(t) < 0 for all time t > 0. On the other hand, x^(t) and X5(t) be- 
ing mass and energy, are both nonnegative quantities for all times t > 0. 
Furthermore, the quantities A, R, and V3 are all positive scalars. Thus, 
the above assertion can not hold. 

The above property of the in-stall reverse flow model is intuitively 
correct, ais the reversal of the flow direction in a real world situation is 
only temporary and can not represent a steady state operating condition. 

In this way when we speak of equilibria in this work we are referring 
to the steady state operation of the compressor (model) in normal unstalled 
operating conditions or in in-stall conditions with normal flow. 
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5 .3 Generation of Equilibria By Continuation 

Let's look at a stationary point of Equation (3.1). This makes the 
left-hand side equal to zero, 

0 = f(x,y) . (5.1) 

Treating y as a parameter to be varied, the vector x can be viewed as the 
dependent variable. Continuation methods deal with the problem of solving 
(5.1) for the state vector x as a function of the parameter y; they do so 
by utilizing results from differential calculus and topology. In general 
y e R m but here let's assume for the sake of simplicity that y e r 1; as be- 
fore assume x e R n . 

Differentiating (5.1) with respect to y yields 

F(x,y>|^ + -|^(x,y) = 0 (5.2) 


where. 


F: 


i! 

3x 


(5.3) 


It is clear that at a particular x as long as F is one-to-one, successive 
values of x can be obtained as y is varied from some initial value yg. 
Typically the initial value of y is chosen at any point which is computa- 
tionally convenient; the desired value of (x,y) is then arrived at by "con- 
tinuation" of solution from (xg,yg). (In the language of Chapter 3, the 
stationary points at which F is invertible are said to be hyperbolic. The 
points at which F becomes singular are called bifurcation points; at these 
points more than one x satisfies (5.2), for a given y. Singular points at 
which the rank deficiency of F is one are referred to as simple bifurcation 
points — otherwise, they are called general bifurcation points. 


- 27 - 



ScientificSystems 


5 .4 BISTAB: Computational Tool for. Bifurcation and Stability Analysis 

In order to trace the equilibrium solutions of (5.1) for the compres- 
sor models (US) and (PS) the computational software package BISTAB was em- 
ployed; the package was developed by Scientific Systems, Inc., specifically 
for this type of analysis. BISTAB contains three user interface FORTRAN V 
functions and a number of supporting functions and subroutines. They: 

(i) Trace connected bifurcating branches of solutions to Equation 
(5.1), starting from some given initial point. 

(ii) Provide stability analysis of the branch curves and report the 
linearized stability type. 

(iii) Spot Hopf, simple, and general types of bifurcation. (At present 
only solutions starting from simple bifurcation points can be 
traced.) 

(iv) Read in the files generated by the bifurcation and stability ana- 
lysis routines, and sort the points into curves of different 

bifurcation or stability type, for plotting. We will not engage 
in any in-depth analysis of the BISTAB package here; instead the 
reader is referred to a recently published description of the 
package by E.F. Wood, J.A. Kerapf, and R.K. Mehra (1984). In the 
following section the equilibrium loci generated by BISTAB are 
presented. 

5 .5 Equilibrium Solutions for (US) and (PS) 

The unstalled system (US) equilibrium solutions are shown as a func- 
tion of the parameter A, the downstream flow area in Figures 5.1 through 

5.9. Each figure depicts a state variable versus the parameter A. The 

symbols O and * on the diagrams stand, respectively, for the stability 
and instability of the system (US). More precisely, the points labeled as 
□ and * correspond to hyperbolically stable and unstable equilibria, 

respectively. 

As will be discussed later, the question of stall must be related to 
the stability of the underlying dynamic system. Thus we seek the stability 
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behavior of (US) as the parameter A varies. We show that the use of A as a 

proxy for system stability may lead to some inconsistencies. We illustrate 

2 

the point by fixing PN at 60% (e.g., Figure 5.1). 

It can be seen from Figure 5.1 that the stability effects of the para- 
meter A on the system (US) is dependent on the location of the state in the 
state space. Thus as the area A decreases, say, from 28 to 25, the stabil- 
ity status of the system changes three times: at points 12, 23, and 34. 
In the regions 1 and 3, the system is stable, and in the regions 2 and 4 it 
is unstable. The points 12, 23, and 34 represent "border" points. 

This means that a decrease in the flow area may or may not be desta- 
bilizing — at least in a local sense. For example at an (unstable) equi- 
librium point in region 3 which is "close" to the point 23, a "small" de- 
crease in the flow area pushes the equilibrium point towards "inside" re- 
gion 3. In this case a decrease in A is stabilizing — or, conversely, an 
increase in A is destabilizing. The above example suggests that a higher 
setting for the flow area A is not always associated with a higher degree 
of stability. However, for (stable) equilibrium states in region 1, an in- 
crease in A is always (further) stabilizing — as Figure 5.1 indicates. 

In the same manner, (the stable) equilibrium points that are located 
in the neighborhood of the point 23 of region 2, are pushed toward stabil- 
ity by a "small" decrease in A; conversely, they are destabilized by a 
small Increase in A. Changes in A have the opposite effect on the stabil- 
ity characteristics of the neighboring equilibria of the point 12 of the 
region 2. However for points in region 4 decreases in A do further desta- 
bilize the system — and vice versa. 

Now, it can be shown that the flow area A and the corrected flow h^ 
move up and down together — in this same region of interest. (In fact, 
the relationship is almost linear.) Thus the above discussion applies 


This discussion was presented in SSI report (June 1 - Aug. 30, 1984) to 
NASA Lewis. 
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equally well with respect to the corrected flow — viewed as a parameter. 
In this way, for example, an increase in the corrected flow is not neces- 
sarily associated with "more" stability and vice versa. Therefore, it 
seems clear that the observed (or computed) value of the corrected flow 
does not — in a one-to-one manner — determine the system's stability rat- 
ing. 

We conclude our discussion of the subject by noting that the above 
reasoning applies equally well to the uncorrected flow &2(= X 3 ) . In Chap- 
ter 7 we show that the situation becomes even more complicated when we en- 
counter non-steady state "off-design" conditions. We will suggest a more 
general approach to the problem which will alleviate the situation. 

Finally, we wish to point out that the in-stall equilibria for this 
range of parameters turn out to be stable hyperbolic. The corresponding 
figures have not been included, for the sake of brevity. (Some equilibria 
are shown in Figures 6.1 through 6.3 in the next chapter.) 
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6. BIFURCATION - CATASTROPHE LOCI 


6.1 Introduction 

In the previous chapter equilibrium loci of the compressor model in 
(US) and (PS) modes were described. This chapter emphasizes the bifur- 
cation loci. It is shown that the unstalled compressor exhibits Hopf 
bifurcation. The resulting family of limit cycle solutions are propagated 
by parameter continuation; both stable and unstable limit cycles are 
observed. Before the unstable family disappears, it further bifurcates 
into 3 families. Frequency, amplitude and centroids of the cycles are 
plotted in various state-parameter spaces. (The analysis is presented in 
Section 6.3.) 

Section 6.2, however, presents bifurcation points of a different type: 
those corresponding to the switching of the models. Strictly speaking, 
such points should perhaps be labeled as "induced bifurcation" points since 
the jumps are induced by switching of the models. Hysteresis is shown to 
exist in the stall/recovery process. 

6 .2 Bifurcation of Solutions and Stall/Recovery 

Let's recall from Section 4.5 that the overall compressor model con- 
sists of three interlocking models (US), (PS), and (RS). The models are 
switched based on the threshold values of the corrected flow: t r (PN) , 
t s (PN), 6 , and e (e.g., Figure 4.1). 

In our analysis we have taken 5 = e = 0, indicating that at nearly zero 
but slightly negative flow the (PS) model is replaced by (RS) ; at zero flow 
(RS) is replaced by (PS). The values of t r (PN) and t s (PN) are given in 
Appendix A. 

Figures 6.1 - 6.3 depict the loci of equilibria and jump points for the 
values of corrected fractional speed of 60%, 70%, and 100%, just to cite 
three examples. Let's note that for example in Figure 6.1(a), we have the 
bifurcation diagrams for the (state, parameter) combinations (x^.A), 
i = 1,5; each x^ represents a state variable. For example, X 3 represents 
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the airflow rate W2 through the compressor. 

Labels (USE) and (PSE) are used for indicating the equilibria of the 

unstalled and in-stall positive flow models, respectively. The dashed 

• 

lines are the bifurcation points, and represent catastrophe points of the 
process using Thom's terminology (1974); the arrows indicate the direction 
of the movement of equilibria and jumps. Thus, for example, on the (x3,A) 
diagram of Figure 6.1(a), as the downstream flow area A of the unstalled 
model is decreased from, say, 32, the flow X3 continuously drops from a 
value of approximately 30 to that of 25 corresponding to A 3 25. At this 
point there is a sodden drop in X3 to the stalled value of X3 3 19. 

This process can also be described, for example, in terms of the 
compressor pressure ratio f} (output 6 of Figure 6.1(b)) and/or temperature 
ratio f2 (output 7 of Figure 6.1(b)). Thus at the stall point 
corresponding to A 3 25 the pressure ratio falls from approximately 1 .63 to 
that of 1.34 and the temperature ratio rises from approximately 1.18 to 
that of 1.20. 

Once the compressor enters the stall mode, recovery requires "extra 

effort" to take place. This is the result of hysteresis in the process. 

As A is now Increased from the stall value of 25 it must reach the value of 

approximately 33 before recovery can take place. At this point the flow X3 

jumps from the stall value on (PSE) of about 27 to the unstalled value of 

about 31 on (USE). Or in terms of pressure and temperature ratios, the 
* 

pressure ratio P3/P2 Jumps from the stall value on (PSE) of about 1.38 to 

* 

that of about 1.59 on (USE); the temperature ratio T3/T2 drops from the 
stall value of about 1.17 to that of 1.15. Figures 6.2 and 6.3 contain 
the same type of information for PN » .7 and PN ■ 1.0. (We did not feel it 
would be useful to include other speed values here — though they were 
computed. ) 

It is noted that the stall/recovery hysteresis loop grows as the 
corrected speed PN is increased. This is due to the fact that the stall 
thresholds increase with PN, while the in-stall characteristic remains 
relatively unaffected; this increases the unstall-to-stall gap at the point 
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of stall. However, the stall-to-unstall gap may or may not Increase — 
depending on which segment of the unstalled characteristic the recovery 
from stall to unstall takes place (e.g., compare gaps for PN = .6, .7, and 
1 on Figures 6.1(a), 6.2(a), and 6.3(a)). 

Figures 6.1 - 6.3 also illustrate that higher PN values require larger 
throttle variations, i.e., the width of the hysteresis loop increases with 
PN. In this way, for compressors with small allowable variations in 
throttle ares, recovery from rotating stall becomes a difficult problem. 
"Small" discontinuity in pressure change at stall is often associated with 
the compressor operating with one or more stall cells that do not cover the 
total height of the annulus. This is known as part span stall. For this 
type of stall the hysteresis effect is relatively small. On the other hand 
"large" pressure drops (associated with large hysteresis effects) charac- 
terize the full span stall. Figures 6.1 and 6.3 illustrate these cases. 

We have been speaking of stall/recovery in the context of bifurcation 
of solutions. Figures 6.1 - 6.3 represent loci of equilibria and jumps. 
Thus the equilibrium loci describe the steady state behavior, not the tran- 
sient. In other words the notion does not involve the so-called 
"off-desig'n" situations. For this reason the stall and recovery process 
described above is perhaps best viewed as local (static) stall and local 
(static) recovery. 

In order to capture the global transient or dynamic nature of the pro- 
cess, a different notion is required, namely, the domains of attraction of 
various equilibria and limit cycles. Though this subject is clearly beyond 
the scope of this Investigation, we do nevertheless touch upon the subject. 
However, this is deferred until the following chapter in which the 
stall/recovery strategies are discussed; in that chapter parameter sen- 
sitivities are also discussed. 

We note that the previously described local and global notions are 
intimately related — as the former provides the core for the latter. More 
specifically, the equilibrium loci provide the attracting and repelling 
poles for the trajectories of the dynamic process; thus they help shape the 
various zones. 
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We close by mentioning that the terms dynamic and static have been used 
in the literature of stall in different contexts. In particular, Greltzer 
(1980) discusses notions of static and dynamic stability. Roughly 
speaking, his description of dynamic stability corresponds to our defini- 
tion of (local) static stability. 

6 .3 Hopf Bifurcation 

Conditions for existence of Hopf bifurcation were given in Section 3.4. 
Here we describe the computational and algorithmic aspects of the problem, 
together with graphical illustrations. 

Let’s recall that we are Interested in (oscillatory) solutions of the 
nonlinear dynamic systems 


dx 

dt 


f(x,u) 


with the property that 

“ *t+T " 0 


(3.1) 


( 6 . 1 ) 


for some xg e RQ and T e [0,«), where we define to be a 

solution to (3.1) starting from xq, with the parameter u set at ug. 

Let’s define y e RO+l as 

y ■ U.ul’ (6.3) 

and Introduce an arc length parameter i, giving 

y<*) - tx(i),p(i)] ' . 

Then (3.1) and (6.1) become 

f(y(*0)) - $(y(*o)) - $ t +T(y<*0> ■ 0 (6.4) 

where f e R n . 

We are Interested in finding n+2 unknowns (T(4g), y(ig) satisfying 
(6.4). Suppose that for some ug we know that there is a limit cycle solu- 
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tion; then we seek the (n+1) vector (x 0 ,T) satisfying n relations in (6.4); 
thus we need a further condition. Chua and Lin (1975) describe a procedure 
for fixing one element of xq to carry out the necessary predictor corrector 
iterations for hoaing in on the cycle. 

The generic software package that we have developed contains a para- 
meter continuation feature that generates a family of limit cycles as a 
function of the continuation parameter p , starting from some easily 
obtained known solution corresponding to some initial parameter value tig.-* 
(We will refer to this generic package as the Hopf package.) Thus starting 
from some initial point "near" a limit cycle, the Hopf routine will close 
in on the exact location and period of the cycle, for the particular value 
of the parameter u that was selected. The continuation part of the routine 
consists of updating the continuation parameter u and finding the 
corresponding new limit cycles that are created. (See Carroll and Mehra 
(1982) for a description of BACTM.) 

Just as in the case of stationary equilibria, a limit cycle may or may 
not be stable. This means that in the case of a stable limit cycle, ini- 
tial points beginning near the limit cycle will converge to it; the con- 
verse is true of unstable limit cycles. Hlrsch and Smale (1974) present 
conditions for stability of limit cycles. The issue is related to the sta- 
bility of the matrix G 


G 


3 

:* . ■ ■- 
3 XQ 


(6.5) 


Here we are Interested in the stability of the discrete mapping $ Q — as we 
wish to compare the values x(tg), x(tg + T) , x(tQ + 2T) , . . . x(tg + nT),... 
for various nonnegative integers n.^ The precise condition for the stabi- 


This is the generic form of BACTM. BACTM is the acronym for Bifurcation 
Analysis and Catastrophe Theory Methodology. The software was developed 
by Scientific Systems, Inc. for the analysis of flight dynamics problems. 

is defined as the composite mapping $.$. ..., repeated n times. 
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lity of a limit cycle Is that all but one eigenvalues of G must be less 
than 1 in magnitude, while the remaining eigenvalue is always 1. 

Unstable limit cycles are not easily observed in nature or in time 
simulation experiments. However, their presence does change the underlying 
dynamics. In the context of our investigation an unstable limit cycle can 
"carry” a trajectory into the in-stall (normal or reverse flow) region, or 
into the unstalled region — depending on the exact position of the cycle 
and the Intervening control actions adopted. 

Returning to our discussion of compressor stall, we are Interested in 
determining whether the compressor exhibits any oscillation of the type 
described by the Hopf theory. This type of oscillation is often asso- 
ciated with surge in compressors (e.g., Greitzer (1982)). The present 
investigation appears to be the first published work indicating that 9uch 
oscillations are due to Hopf bifurcation of solutions. 

The unstalled compressor model (US) exhibits Hopf bifurcations at cer- 
tain critical values of parameters. We have observed this behavior for 
every given corrected speed PN value (of course at different parameter 
settings). We have verified the oscillatory behavior by time simulations. 
Figures 6.4(a) through 6.4(f) show the development of a stable limit cycle 
for the corrected speed PN » 0.6. It can be seen that before the limit 
cycle disappears it becomes unstable. The stable limit cycle family 
reaches its peak amplitude for A = 27.2. Figure 6.4(f) indicates the 
period of the limit cycles as the parameter A varies. The peak period of 
the limit cycles is reached for A a 27.32; the corresponding period Tp^^ = 
16.6 x 10~3 seconds. The minimum period of about 16.36 x 10“^ seconds is 
reached for A ■ -26.9. 

As A is increased from the value of approximately 27.2, the stable 
limit cycle amplitudes decrease. For A = 27 .3892 and beyond the stable 
limit cycles disappear and a family of unstable limit cycles is born. 
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s: stable Halt cycle 

u: unstable Halt cycle 
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Figure 6. 4(c), (d) . Development and Disappearance of Stable Limit Cycle 

Projected on (x^, A) Plane, i » 3,4. 

(State 3 ■ u^, State 4 - W^) 
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Figure 6.4(e). Development and Disappearance of Stable 
Limit Cycle Projected In (x 5 ,A) Plane. 
(State 5 - (WT) 3 ) 
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Figure 6.4(f). Period of Limit Cycles Versus A. 
(State 6 • T ■ period) 
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6.4 Global Bifurcation of Limit Cycles 

Equation (6.1) of the previous section describes the property that 
limit cycle solutions must possess, namely, the existence of a pair (xq.T) 
such that a trajectory of (3.1) satisfies 

♦t(xO»Uo) “ ♦t+T( x 0>l 1 0) " 0 (6.1) 

for all times t > 0. Equation (6.5) was used to describe the local stabi- 
lity of a limit cycle ‘’starting" at xq and having the period T. Just as in 
the case of stationary equilibria, at some critical parameter value, say 
uq» Equation (3.1) may admit a multitude of solutions, say m of them, with 
the property 


*J (x 0’ u c> " *{+T 1 (x 0’ w c> " 0 < 6 * 6 > 

for "initial condition" x* and period T^, i ■ 1 m. 

Each limit cycle would then have its local stability property as 
described by the corresponding stability matrix G*, i - 1 


Gi 


3 (frt(*0»V c ) 

3 xj 


(6.7) 


If m represents the total number of limit cycles, then their determination 
represents a global solution to the problem. The situation is depicted in 
Figure 6.5. In this figure a total number of 3 limit cycles in the solu- 
tion space X have been depicted for the parameter value u c . 

Figures 6.6(a) - (b) depict the development and disappearance of a 
family of stable limit cycles — giving rise to an unstable family which 
subsequently bifurcates into three families of limit cycles; the con- 
tinuation parameter is A, as before. Only 4 of 6 state variables are shown 
— for the sake of brevity. Vertical bars indicate the amplitude of the 
limit cycles. As can be seen from Figures (a),(b),(c), and (d) as Che 
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parameter A is increased from its initial value of about 27.3890 the ampli- 
tude and the period of the stable limit cycles continues to shrink. 
Approximately at A • 27.3892 and beyond the limit cycle becomes unstable 
while decreasing in amplitude. Further, increase in A to a value of 
approximately 27.3895 causes an increase in the number of limit cycles 
(from 1) to 3. As can be seen from the Figures, this condition persists 
until A exceeds a value of approximately 27.38955. At that point and 
beyond once again a unique (unstable) limit cycle family is generated. It 
can be seen from Figure 6 . 6 (d) that the three limit cycles not only have 
different amplitudes and "locations" but that their periods are also dif- 
ferent. We note that the above graduations in A are too small to be experi- 
mentally verifiable, even though the software was calibrated for a 10 ~^ pre- 
cision level. The point of the above discussion is to point out the theore- 
tical possibility of obtaining 'multiple surge cycles' for some critical 
parameter values. In fact ve do not claim to have discovered all such cy- 
cles. A further investigation is warranted to either validate or . refute 
this possibility in a number of realistic models. 

Ue close by noting that the above solution classes represent a fairly 
complex dynamic phenomenon, Involving stable and unstable limit cycles of 
varying quantity, amplitude, period, and location — all derived as a func- 
tion of the parameter A. In this way it can be observed that the usual 
one-to-one association of the values A with "stability" or "instability" is 
at best an ambiguous notion. This same observation can be expected to 
extend to other parameters. 

6.5 Time Simulations 

Here we will present several simulations of the process for various 
values of the parameters PN and A. This will also serve as a confirmation 
of some of the results obtained in previous sections by other means. 


and output 


Figures 6.7(a) , (b) , (c) , and (e) present the time profile of three 

dW« p* 

variables: state variable X 3 :- —2. , output variable Fit- 3 . 

T, dt ?2 

variable Fo ■- _ . Other variables have similar properties and have not 

2- T7" 

been sketched, for the sake of brevity. 
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Figure 6.7(d). Phase Portraits 
of the Stable Limit Cycle in 
Figure 6.7(c). Cycle is Projec- 
ted on (xi,H|) Planes for i * 1,5 
Hi * Corrected Flow Rate Through 
Compressor. 
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Figure 6.7(a) presents a trajectory with growing oscillations which 
begins in (US) and ends in rotating stall (PS). An increase in A to the 
value of 26.7 restores equilibrium in the unstalled region; this is 3hown 
in Figure 6.7(b). A further increase in A, however, brings about oscilla- 
tory behavior. Comparing this case with that illustrated in Section 6.5 
(e.g.. Figures 6.4(a) through 6.4(f)), we see that these oscillations are 
indeed limit cycles produced by Hopf bifurcation of equilibrium solutions. 
It can be observed from the plots that both amplitude and frequency of 
oscillations agree very closely. 5 A Phase Plot of the time function of 
Figure 6.7(c) is shown in Figure 6.7(d); it shows how a trajectory 
approaches the (stable) limit cycle. A still further increase in A to the 
value of approximately 27.5 further restores the unstalled equilibrium. 

Thus from the point of view of recoverability it makes a difference by 
how ouch a certain nominal value of A is changed and in what direction. 
Both decreases and increases in appropriate ranges can remove surge type 
oscillatory behavior. It should be observed that the implications for 
surge/recovery due to area changes described in the previous paragraphs (or 
changes in any other parameters) are dependent on the particular trajectory 
in question, with respect to its location in the trajectory space. Thus 
further characterization and analytical machinery is required. In the next 
chapter we will address this problem in part. Finally, surges in the state 
variables rate for corrected speeds PN of .7, .8, .9 and 1 are listed in 
Figures E(a) through E(d) in Appendix E. 


^It is clear that not every oscillatory motion corresponds to a limit 
cycle . 
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7. DOMAINS OP ATTRACTION FOR OFF-DESIGN STABILITY, STALL/ RECOVERY , AND 
SENSITIVITY ANALYSIS 


7.1 Introduction 


The problem of compressor stall is ultimately related to the stability 
behavior of the compression system under conditions other than the normal 
steady-state operating environment. In Section 5 we discussed the loci of 
solution equilibria; and the stability issues that were discussed referred 
to the local properties of equilibria. However, a given trajectory may or 
may not be 'near* any equilibria; its stability behavior depends on its 
overall position, in a global sense. 

In order to obtain a complete classification of trajectories we need 
to compute the domains of attraction of equilibria. This problem, however, 
is beyond the scope of this investigation. But in order to highlight its 
usefulness for stall/recovery we discuss the issue in Section 7.3. In 
Section 7.5 we show the applicability of the domain of attraction for use 
as a criterion of stall. This appears to offer a natural generalization of 
the notion that the (low) values of corrected flow can be used as proxy for 
entering stall, e.g.. Section 7 .A . Partial computations of domains of 
attraction are presented in Section 7.6. However, in Section 7.2 we first 
define three types of trajectories according to the 'severity' of their 
s tall/ recovery, and give some examples. 

7 .2 First Stage, Second Stage, and Third Stage Stall/Recovery 

We have .seen that the total dynamic behavior of the compressor in 
various conditions can be modeled as three Interlocking dynamic models 
(US), (PS), and (RS). Let's start with a compressor in normal steady state 
unstalled condition and perturb the steady state operating point, say xq, 
to that of * xq + 5xq to create an off-design condition. A trajectory 
based at xj and under the control regime U(*) is defined as $ t (xi,u). 
Depending on the stability properties of xg and the amount of perturbation, 
the resulting trajectory may (i) only involve (US), or (ii) involve (US) 
and (PS), or (Hi) involve (US), (PS), and (RS). 
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One can Chink of off-design transient conditions which, using the 

available control action u, can be cleared without having the compressor 

enter into either of the two stall conditions as being somewhat benign. We 

$ 

define such trajectories as being first stage recoverable. On the other 
hand, if despite all the available control action the trajectory still 
enters stall but can be cleared without entering reverse stall, such a 
stalling trajectory is defined to be second stage recoverable. Finally, if 
the stalling trajectory can be cleared only upon entering (and leaving) 
both stall conditions, the stalling trajectory is defined as third stage 
recoverable . 

The above terminology can be justified by noting that (i) physically, 
every time a compressor stalls it stands a certain chance of being damaged, 
perhaps due to some unmodeled factors and, (li) from a modeling point of 
view a first stage recoverable trajectory will Involve only one dynamic 
model, thus involving the least amount of dynamic complexity, etc. 

Figure 7.1(a) shows a trajectory that begins in the unstalled region 
and eventually settles in a rotating stall equilibrium. This stalling 
trajectory can be cleared by a combination of a decrease in PN (to 60%) and 
an Increase in A (to 32). Figure 7.1(b) shows that the trajectory is first 
stage recoverable. On the other hand, Figure 7.2(a) shows a different 

stalling trajectory. It turns out that for this trajectory no combination 
of parameter setting exists that can keep it from entering stall. Figure 

g 

7.2(b) shows that this trajectory is second stage recoverable. 


^These terms are defined intuitively only; one can refine these definitions 
by introducing a mathematical notation. 

For the sake of brevity only x, component of the trajectory is shown. 

8 J 
Here we limit ourselves to fixed parameter settings. 
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(a): PN-.7, A=30.25, x 3 «))-30.0 (b) : PN-.6, A-32.0, x 3 <0)-30.0 


Figure 7.1. Example of first-stage recoverable trajectory. 

Using combination of PN and A, a stalling tra- 
jectory (a) is recovered without entering stall (b). 



Figure 7.2. Example of second-stage recoverable trajectory. 

No combination of PN and A exists which can keep 
trajectory from entering stall before recovery. 
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7.3 Domains of Attraction 


Intuitively an attracting set A for a point x is that which 'attracts' 

the trajectory $ t (x,uo) baaed at x as time increases. 9 Accordingly, the 

domain of attraction of A is made up of all those xs which will eventually 

end up (under the action of the system) at A. A little more notation makes 

this more precise. A closed invariant set A C R° is called an attracting set 

if there is some neighborhood V of A such that ^t( z ) ♦ A as tx», for all 

xeV. The set D:- D $ t (V) is said to be the domain of attraction of A. 
t<0 

(The set A may be just one (equilibrium) point.) 

These simple ideas can be used to advantage for the analysis of stall 
and recovery. Let a normal operating condition of the compressor be 
denoted by xq, corresponding to the nominal parameter no* The elements of 
iiq may be some design parameters (such as the downstream volume V 3 ), and 
some adjustable parameters (such as the throttle area or the rotor speed). 



It is observed from Figure 7.3 that trajectories based at points outside 
Dyg will stall. (For the sake of argument, assume that all those points 


9 


Assume that the parameter u is fixed at 


UQ« 
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not converging to xq will end up in rotating stall.) However, by manipu- 
lating the parameter y, say resetting it to y^, the shape of will 

change to that of Dy^. It is observed that in this case an otherwise 

stalling trajectory such as one based at y will now approach the unstalled 
stationary operating point y. In this way be varying the parameter y over 
the entire feasible parameter space U one obtains the overall recoverabi- 
lity domain as 


Qj i m U Dy (7.1) 

u yeu 

In this way represents the set of all trajectories that can be 

recovered, using the feasible parameter set u. 

The set can be computed and stored off-line, once for a given dyna- 
mic system. Now if we can pinpoint the location of the state x, the 
question of whether the trajectory based at x can be recovered reduces to a 
table look up probelm: determine if the trajectory belongs to Q>\ if so, 

the trajectory is recoverable - otherwise it is not. 

We can utilize the definition of recoverability types (e.g., 1 st , 2 nd , 

and 3 rd stage recoverability) to obtain three types of recoverability 
12 3 

domains : ^>y , ^y , and ^y . 

7 .4 Transient Analysis and Corrected Flow 

(a) In this section we first note that in non-transient situations the 
uncorrected flow rate through the unstalled compressor is a good proxy 
for its state. In fact in steady state, the uncorrected flow rate X 3 
by Itself determines the remaining state variables and the corrected 
flow* This will in turn determine various pressure and temperature 
ratios* This can be seen by solving the steady-state state vector in 
(US) 


0 - F(x,y ) (US) 

where F represents the RHS of Equations (U1)-(U5) in Chapter 4. 
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The particular form of functional relationship is somewhat complicated 
as it involves a series of substitutions involving nonlinear 
relationships. 

On the other hand, as far as stability is concerned, we have seen that 
there is no one-to-one (monotone) relationship between the value of 
the corrected flow h^ (or the flow area A) and the local stability of 
equilibria. Thus even for local stability determinations the 
corrected flow value does not contain the necessary information. The 
reason is that the corrected flow variable h} is an output variable of 
the system (US) the value of which is determined in a unique casual 
fashion, by the state vector x. 

The question of whether a trajectory is in-stall must therefore be 
related to the location of the state vector x in the state space 
rather than that of the output level(s) in the output space. It is 
clear from Equation (4.10) that there are many state vectors x that 
will give rise to the same h}. This is depicted in Figure 7.4. 


l®This was done in SSI's Progress Report (Jan. 1 through May 1, 1984) to 
NASA Lewis. 
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sec of admissible x set of admissible h} 

Figure 7.4 The corrected flow hi does not uniquely determine the 
8 cate x. 

Therefore the output value hi does not give complete Information about 
the state x. 

(b) It can be noted (for example, from the analysis of Section 7.6) that 
the situation Is even more complicated in case of off-deslgn 
(non-steady state) conditions. In particular, no one state variable 
can be used to determine the remaining state variables; knowledge or 
estimation of individual state variables Is required. Furthermore, 
the exclusive use of corrected flow rate (or A) for stability rating 
in this transient case is even less acceptable. 


7 .5 Domain of Attraction as Stability Set and Criterion for Stall 

The definition of stall used in this work** is that spelled out in 
Section 4.6. In particular, for a given corrected speed PN there is a 
threshold stall value, t s (PN) for the corrected flow. As the trajectory of 
the system evolves over time so does the corrected flow h}. If at any 
point in time hi falls below t 3 the unstalled compressor model (US) is 
replaced by (PS) , and whenever hi equals or exceeds the recovery threshold 


l^This definition is the same as that spelled out in the Statement of Work 
(e.g. Wenzel, Bruton (1982)). 
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t r (PN) the compressor becomes unstalled. Thus the models are switched 
based on the scalar output variable h} - even though, as we have seen, this 
variable does not contain the entire stability information. In this way 
this choice of stall criterion may for example lead to stalling of some 
trajectories that would otherwise recover, or in some cases may provide a 
'late' signal of oncoming stall. 

Based on the earlier discussions, it appears that a more general and 
reliable criterion of stall would be provided by the domain of attraction 
of some 'normal' operating point (or operating region.) 

Let xg denote some nominal steady state vector corresponding to the 
corrected spped PNg, and flow area Ag. Let D(PN 0 , Ag) denote the 
corresponding domain of attraction for xg. Then one way to decide how to 
switch the models is as follows : 

(I) If x(t) e D( PNq , Ag) , the trajectory will recover, approaching xg 
as time Increases. 

(ii) If x(t) t D( PNg , Ag), then replace (US) by (PS) provided h^O. 

(Hi) while using (PS) at any time hi<0, switch to (RS) until the 
flow becomes nonnegative again at which time switch to (PS). 

(iv) If while using (PS) the trajectory enters D(PNg, Ag), then switch 
to (US). 

Figure 7.5 illustrates the procedure. 
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x £ D(PN 0 ,A 0 ) 



The previous rule a sy be viewed as being applicable to passive 
stall/recovery. That ls t In this case no control action Is taken when the 
trajectory goes outside D(PNo, Aq>* As we illustrated previously, one aay 
nanlpulate the control parameters to obtain D, so u to atteapt to recover 
a particular trajectory. Then the switching rule may be based on this new 
set or some variation of it. We will not elaborate on this point here. 
Furthermore, we have not made any distinctions about different stages of 
recovery here although this can be easily Incorporated. 

In closing we note that the set D(PNq, Aq) Is contained In R n . This 
la in contrast to the scalar output t 8 (PN) rule of Wenzel, Bruton (1980). 
It appears that the use of the set D(PN) in place of the scalar threshold 
value t a (PN) represents a natural generalization of the latter approach. 

7 .6 Some Experimental Results for Stall/Recovery and Sensitivity Analysis 

Some partial computations of the domain of attraction were carried out 
to illustrate the usefulness of the approach. The five dimensional space 
was projected on the ( 23 , X 5 ) plane. The state variables xj, 13 , X 4 were 
fixed at values 1.8, 940, and 0.2, respectively. These values are dif- 
ferent from their corresponding equilibrium values; they correspond to some 
arbitrarily chosen 'off-design' condition. This aay represent some per- 
turbed values of some normal operating condition. The nominal area A was 
fixed at 30. 25. 


- 71 - 



Scientific Systems 


The computations were carried out for two sets of parameters: (1) 

nominal, and (ii) perturbed (increased from nominal value by 10Z in most 
cases)* The particular parameter no m i na l values chosen and their perturbed 
values are: 

nominal parameter values 

perturbed parameter values 
(For definition of parameters see Appendix A and Wenzel, Bruton (1982).) 

Figure 7.6 shows the domains of attraction for the unstalled system 
equilibrium (USE) and stalled system equilibrium (SE) 12 . In each diagram 
USE1 and SE1 refer to the two equilibrium values for the nominal case 
(shown in solid lines), while USE2 and SE2 refer to the perturbed case 
(shown in dotted lines.) The arrows on the figures point to the stall 
region. It is seen from Figure 7.6 that the domains of attraction of the 
two equilibria split the (X3, xj) plane into several regions. For example 
Figure 7.6(b) shows the stall and recovery regions for the zero flow effi- 
ciency values no ■ 0.1 and no ■ 0.05. Focusing on the nominal case 

no * 0.1 (solid separatrlx lines), let's fix attention on the 'horizontal' 
line: X5 * 60.0 and X3 free. We see that for the flow rate X3 between 

about 5 and 25 the resulting trajectories will all end up in stall. Both 
higher and lower values of X3 with respect to this range will result in 
recovering trajectories. We can also see that by fixing X3 ■ 30.0 and 

allowing X5 to vary, five different regions can be very complicated. We 
may have various 'pockets' of stall scattered around the recovery regions - 
a situation which may not be discernable from equilibrium or local stabi- 
lity analysis. 


A - 30.25 
R - 639.6 
L - 9.64 x 10 -4 
no ■ 0.1 

A - 27.5, and 27. 3j 
R - 703.56 
L - 1.0604 x 10“ 3 
no * 0.05 


* 2 Here we make no distinction between various types of stall and recovery 
discussed in the previous section; only the final point of convergence is 
considered. 
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Now by varying the parameters, the shapes and sizes of various regions 
change - providing information on which a set of trajectories may be reco- 
vered by changing the parameters in their feasible sets* For example by 
changing the zero flow efficiency variable hq from its value of 0.1 to that 
of 0.05 we obtain the new separatrix shown by the dotted lines. In this 
case we note that the parameter change increases the stall region, but 
slightly. This provides a method of sensitivity analysis for the impact of 
various parameters on the stall/recovery regions. Due to the fact that, in 
general, as a result of parameter changes we obtain a complicated pattern 
of change in the various regions, sensitivity ranking of the parameters is 
at best an ambiguous notion. Wenzel and Bruton (1982) provides a sen- 
sitivity analysis of some parameters but their method of analysis uses the 
turbine area A as the stability ranking mechanism. However, as can be seen 
from Figure 7.6(d), the value of the area A does not seem to provide a 
clear cut measure of recoverability - especially in some critical ranges. 

Figures 7.6(a), and (c) provide the sensitivity of the stall/recovery 
regions to the changes in the parameters \ and R. They have both been com- 
puted based on a 10Z Increase in the nominal parameter values. Of the two, 
it appears that Increases in R have a greater Impact on increasing the 
stall region. 

In closing we note that the above analysis is not intended to provide 
specific recommendations for stall and recovery; rather, our goal has been 
to highlight the potential usefulness of the techniques for providing spe- 
cific strategies for stall avoidance and recovery. 
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8. CONCLUDING REMARKS AND SUGGESTIONS FOR FURTHER RESEARCH 

This Investigation illustrates that analytical methods employed in 
nonlinear control theory are quite appropriate for the analysis of 
compressor stall phenomena* In particular the methods of bifurcation and 
catastrophe theory when combined with a state variable formulation can help 
formulate a general purpose methodology which provides insight into the 

stall problem. The approach can provide an across-the-board understanding 

✓ 

of various models - all within the same analytical framework. 

Though we have presented some sample computations of domains of 
attraction of equilibria, this objective was beyond the scope of this 
investigation and therefore remained Incomplete. As the discussion of 
Chapter 7 indicates, the accomplishment of this objective will provide 
valuable Insight into questions of: parameter sensitivity analysis, 

regions of stall and recoverability from stall, types and availability of 
feasible control action, just to mention a few. 

Switching of the models in this report was based on the threshold 
stall values t s (PN) of the corrected flow rates. As we argued in Chapters 
6 and 7 this criterion may lead to incorrect conclusions. As a result in 
Chapter 7 we proposed a new procedure for switching of the models - 
involving the use of domains of attractions. Further investigation into 
this approach should be carried out. In particular, a careful matching of 
the domains of attraction of the stalled and unstalled models must be made. 
Otherwise a constant back and forth switching of the models in some regions 
may result. On the other hand if the models are treated as given, then 
these potential collisions of the attraction domains may result in a new 
type of 'surge' phenomena. 

The development of state space models invariably brings about the 
question: are the state variables directly measurable quantities? 

Typically state variables represent important 'internal' variables, not all 
of which are measurable. Other related variables (e.g. output variables), 
usually smaller in number than the state variables, may be directly and 
accurately measurable. In our case, various temperature and pressure 
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ratios and/or absolute quantities may be good candidates for measurement. 
Thus observer-estimator type approaches must be adopted to provide closed 
loop estimation of the state vector. (It would be desirable to construct a 
special type of observer for this problem: that which would be 'sensitive' 
to the boundary of the stall region. ) 

The last but by no means the least recommendation is that modern 
control theoretic identification techniques should be brought to bear on 
(i) identification of in-stall characteristics, and (ii) identification of 
model parameters in the face of noisy data. 
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APPENDIX A 

(Provided by NASA Lewis) 


Pi - 14.8 

(Al) 

T x - 519 

(A2) 

• P1-P2 


Wl ■ lo 

(A3) 


t 

w 2 ■ / (»l - £fc)dt + »2,i (A4) 

o 

(WT)2 - Y / (wj - W 2 (pj) dt + (WT) 2fi (A5) 

r(wt) 2 6 

p 2 , T 2 " X 3(1 “T7972) (A6) 

(WT) 2 

^2 “ y 2 (A7) 



(A8) 



(A9) 


N//9 2 

16042 (A10) 


A-l 
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$ - 7.972 x 10-3 


, W2^ 9 2/ 6 2 


PN 


(All) 


+ -© * 2 + 0.33 


(A12) 


UNSTALLED 


STALLED 


- P 2 


■ 


P 3 " + * 7746 * Pl,2) (A13) 


* /W2/9 2 \ 

r 3 - T 2 * F 2 \i — p y 


T, - T 2 


// P *\.286 

I [*2/ 

1 * \~ ' 


- 1 


(A14) 


w sx » f 3 (pn) 


n - no + F^PN.noH 


(A15) 

(A16) 


W 2 J 9 2 

STALL IP ■■ ■ ■■ - W ST 


W 2 /92 
>2 


RECOVER IF " > 1.1 W SX 


W 3 - / (W 2 
o 


- W 3 )dt + W 2> i 


(A17 ) 


* lr * 

w 2 - l / CP* - P 3 )dt + W 2>1 


(A18) 
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(WT) 3 M / {fj} - W 3 T 3 )dt + (WT) 3>; 


W 3 =* A 

S13 


R(WT) 3 

P 3“ — 


*3 


(WT) 3 

w 3 


A ■ Aj for t ti and t t 3 

- A a - a(t - tj) for C| < C < t2 

* A a - a(t - t^) + 2a(t - t 3 ) for C 2 < t} < t 3 


STALL THRESHHOLD VALUES t s (PN) 


PN 

tg(PN) 

0.6 

24.8 

0.7 

31.2 

0.8 

42.3 

0.85 

47.1 

0.90 

51.2 

0.95 

55.3 

1.00 

59.5 


t r (PN) =* 1.1 tg(PN) 


A-3 


(A19) 

(A20) 

(A21) 

(A22) 

(A23) 
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A 

area 

in 2 

m 2 

dt 

time increment 

sec 

sec 

fi 

performance function (i"l,4) 

- 

- 

KN 

constant describing in-stall compressor 




performance at negative flow 

- 

- 

KP 

constant describing in-stall compressor 




performance at positive flow 

- 

- 

L 

fluid inertance, length/(area g) 

sec 2 / in 2 

sec 2 /m 2 

N 

rotational speed 

rev/min 

rev/min 

P 

static pressure 

lbf/in 2 

N/m 2 

P 

total pressure 

lbf/in 2 

N/m 2 

PN 

fractional corrected compressor 




speed 

- 

- 

R 

flow resistance, A (pressure drop)/ 




A (flow) 

lbf sec/lbm in 2 

N sec/KG m : 

T 

temperature 

°R 

K 

t 

time 

sec 

sec 

V 

volume 

ln^ 

m^ 

w 

weight 

lbm 

kg 

• 

w 

weight flow 

lmb/ sec 

kg/ sec 

a 

constant 

- 

- 

Y 

ratio of specific heats 

- 

- 

6 

rati of total pressure to standard 




atmospheric pressure 

- 

- 

A 

Incremental change 

- 

- 

0 

efficiency 

- 

- 

0 

ratio of total temperature to 




standard atmospheric temperature 

- 

- 

♦ 

flow parameter 

- 

— • 


pressure parameter 

- 

- 


Subscripts : 

0 ambient 

1 fan inlet 

2 compressor 

3 compressor discharge 

i initial 

* denotes quasi-steady parameter at compression component discharge 


The terms within braces { } in the equations are switched when flow reverses as 
follows: If flow is positive, in the normal sense, the upper term within the 

braces are used. 
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STATE SPACE REPRESENTATION FOR STALLED MODEL (PS) 


The state vector x that we have chosen is 


w 2 

Ofr ) 2 

li 2 
w 3 

L(wt) 3 J 


We have from Appendix A 
* 1“*2 


Wi 


R r 


(A3) 


W 2 m S (»i - <S 2 ) dr> W 2 < (A4) 

o 

Upon differentiation and substitution, equations (A3) - (A4) yield 
‘ ^ (WT) > - 


or 


Pl R 

*1 ’ *2 - *3 


Given 


(WT) 2 - " “2 {J?}) dT + ^ WT )2,i 


(Bl) 


(A5) 


we obtain 


( 51 ) 


B-l 
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But 


<WT>2 x 2 
2 W 2 


(B2) 

and 



R(WT) 2 r 
p 2 • v 2 • V 2 *2 


(B3) 

yields 



*2 ■ T QO ' nk ‘ 2 1*2/0 

- 3 ^}) ■ 

<B4) 

Now 



t 

W 3 ■ / (u >2 “ W 3 ) dr + ^3,i 


(A17) 

. P3 

013 - A 

^3 


(A20) 

with 



R(WT ) 3 r 
p 3 ■ V 3 ■ V 3 *5 


(B3 1 ) 

(«T ) 3 X 5 

3 w 3 X 4 


(B2 f ) 


yield 
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. AR j 

x 4 » x 3 - — /xax 5 

(B5) 

Equation 


1 C 

<*>2 ■ T / (Pi - P3) dt + o>2 i 
0 

(A18) 

upon differentiation and substitution obtains 


* 3 " r (*3 iv7 * 5 ) 

(B6) 

At this point we shall concentrate on the 
positive flow direction. Thus 

stalled dynamics in the 

f 3 * 0.837 '<* + - 7746 * ™ 2 > 

(A13) 

where 


♦ - KP <p 2 + 0.33 

(A12) 

_ i2'(9/62 
♦ » 7.972 x 10-3 ■ / pN . * 

(All) 

Using (A- 11), (B2) , (B3), and 


T2 

8 , a — i 
2 519 

(A8) 

* P 2 
TTJ 

(A9) 

we obtain 
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81 *3 

' P1) /5l9^I 

where 

A 7.972 x 10-3 x 14.7 *2 

81 r 

Substituting (B7) in (A- 12) yields 
/8l\ 2 x 3 2 

♦ - KP (~) -rr— + 0.33 

VPN/ 519 xjx 2 

Thus 

2 

2 x 3 

P* =■ g 2 x 2 + g 3 g! KP ~ 

82 ■ « (PW - oTahli * °' 305 % ™ 2 

83 ■ 0-925 T j|-^ 

Finally (B6) becomes 

/l . 2 _ x 3 R \ 

*3 “ 82 x 2 + 8381 ^ ~ x 5 j • 

Equations (A- 20), (B2)', (B3)', and 
t 

(WT) 3 - T / (» 2 T* - <i 3 T 3 ) dr + (WT) 3 a 
o ■* ’ 

yield 

*5 * Y^x 3 T^ * x 5 / x 5 /x 4 ) . 

Now 


(B7) 


(B8) 

(B9) 

(BIO) 

(BU) 

(B12) 

(B13) 

(A19) 

(B14) 


B-4 



ScientificSystems 


is - i 2 (i + (^);~ - 1 - ) 


(A14) 


together with Equations (A13) and (A22) obtain 


(P / (sixi)* 

! , L 1,19 * 0,925 (” 519xix2 + °- 


* x 2 
T3 - — 

-» X! 


Furthermore 

n - no + F4 (Pn.no) ♦ • 

Using (B7), (A15), and (A28), 


- - 33 P»2 )]- 286 


-)• 


(A28) 


(A15) 


* x 2 
T 3 ■ — 
J xi 


/ [l.l9 + 0.925 (kP 5 1 9x]x 2 + °* 33 Pn2 )] “ 1 \ 

V 81 x 3 / * 

' no + F4(PN,no)^j — ' 

* *■" / c t ft 


/519 xjx 2 


Upon substitution for T 3 in (B14) 


x 5 


X 2 X 3 


/ fl .19 + 0.1 
1+1 

v no + 


925 (KP + 0 - 33 PN 2 )]* 286 


F4(PN,no>^j 


x 3 


-) 


/519 X]X 2 


- ^ x 5 /x 5 /x 4 . 


(B16) 


deriving the last state equation. Specializing Equation (B4) to the 
positive flow case we collect the equations for easy reference. 
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Summary of State Space Systems for Positive Flow Stalled Dynamics (PS) 


. p i a 

xi " i£ " m * 2 " X3 


(PSD 


x 2 - Y 


(Ik 

\Ro 


Tl ~ RqV 2 Tl X2 ~ X2 Xi 


X 1 / 


(PS2) 


x 3 


2 

x 3 


1 Z ^ *3 ^ 1 /n „v 1 R 

L «3 81 ® ^ + l g 2<PN) x 2 - r — x 5 


(PS3) 


R 


x 4 - X 3 - A “ /X 5 X 4 


(PS4) 


x 5 - y 


X 2 X 3 




0.286 


1.19 + 0.925 /KP ^81 x 3 -) 2 + 0.33 PN^ 


519 x^x 2 


1 ] 


- 1 


no + ?4(PN,no)pjj 


81 x 3 


/519 xjx 2 



where 


XI 


w 2 

x 2 


(WT ) 2 

x 3 

- 

w 2 

x 4 


W 3 

x 5 


(WT ) 3 


(PS5) 
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APPENDIX C 

STATE SPACE REPRESENTATION FOR THE UNSTALLED MODEL (US) 


Equations defined In Appendix A will be referred to as A1 , A2 
etc. Any new notation introduced shall be defined accordingly. 


We have from Appendix A 


wi 


p l~ p 2 

R o 


t • • 

W 2 ■ / (wi - dr + W 2 i 
o ’ 

(A3)-(A4) yield 

• P 1 R 

X 1 “ I; “ x 2 > *3 

Equation (B4), specialized to the positive flow region yields: 


/ (IL- 

R \ 

T, - X3X 0 

\ Uo 

RoV 2 X 2 j 

1 X 1 / 


Equation (A18) obtains the third dynamic relation in the form 
*3 =* Z (p 3 * *5) 

Here P 3 is obtained from the unstalled case of Equation (A13): 


* 

p 3 - p 2 p l( 


• __ 

ui 2 ^ 62 

6 2 


PN) 


(J 


C-l 


(A3) 

(A4) 

(Cl) 

(C2) 

(C3) 

^-13) 
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Let 

• __ 

A o»2^ 9 2 



denote the corrected flow. 


(C 4) 


We have, from Appendix A and definition of x: 
*2 

02 * 519 X! 

62 “ 14.7 V 2 * 2 


(C5) 

(C 6 ) 


yielding 


h l 


x 3 

84 -== 
/x l x 2 


where 

A 

84 * 


14.7 V 2 
R /5IT" 


(C7-) 


(C8) 


(the constants gj , g 2 and g 3 were defined in previous appendixes. 


Equation (C3) thus becomes 



The fourth dynamic relation is the same as in Appendix R. 


(C9) 


x 4 



/x 4 x 5 


Finally,, the fifth relation is derived from (A19) and 
t . * 

(WD 3 ■ y f (wo T 3 - 0)3 T 3 )dr + (WD 3 ^ 
o ' * 

* / U2^2 

t 3 - t 2 f 2 \ — 



(CIO) 


(Cll) 


(C12) 
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Substituting for the corrected flow and temperature T 2 in terms of the 


state variables yields: 



<K4 , X3 = , PN) 


Thus 

* . 

X5 - y(*3 T3 - 013 T3) 


(C 13) 


(C14) 


- T(x 3 -- F 2 (g 4 


^3 

/ *r* i 


PN) 



^5 

X4 


(C15) 


where 

“3 t 3 



CCI6) 
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Summary of State Space System for (Installed Dynamics (US) 


X 1 


£i 

&o 


R 

M2 


x 2 - X 3 


*2 


*3 


Y 

R 

L 



x 2 /xi 

* V 

V 3. 


x 4 • *3 " y7 /x 4*5 



( 01 ) 

(U2) 

(U3) 

(U4) 

(05) 


where the state vector is x 


X 1 


y 2 " 

x 2 

A 

(vrr) 2 

x 3 


u>2 

X4 


W 3 

L*5 J 

■ 

iiWT ) 3 J 
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APPENDIX D 

STALLED DYNAMICS (REVERSE FLOW) 


Since the procedure is Che sane as ChaC of Che unscalled case or 
Che scalled (posicive flow) model, we will derive Che five sCaCe 
equations wlchouc much explanation. 


As before. 


*1 m s~ - 


Pi R 

Ro ' RqV 2 *2 " *3 


buC 


*2 * y 


t 2 - t 3 


(1± _ R *2\ 

\ R o X 1 ~ ^2 *1 / 

(l - - t.. \ 

0.7972J 


x 3 T 2 


T3 


(WD3 X5 
" *4 


_ (D2 /02 

♦ - 7 - 572 !°' 3 if-pr 


(1) 2 / 0 2 


(Bl) ,(C1) 

(D2) 

(D3) 

(D4) 

(D5) 


In-stall 


calculation: 
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0 2 


T3 

Tl 


6 2 


£3 

?! 


P3 



x 5 


T3 


£5 

*4 


Equations (D6)-(D10) yield 


hr 


where * 


<i»2 / 0 2 


86 


x 3 

/ X4X5 


86 


P1V3 

r/T7 


Substitution from (Dll) into (D5) yields 


87 *3 

*P m pjj . 

rn /X4X5 

where 

g 7 - 7.972 • 10“3 g6 
Thus (D3) yields 


T 2 


£5 

*4 



88 

PN 


/X4X5 / 


(D7) 

(D8) 

(D9) 

(DIO) 

(Dll) 

(D12) 

(D13) 

(D14) 

(D15) 


* Subscripts 6,7,... in g's are the continuations of those used in 
previous appendixes. 
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87 

88 “ 0.7972 


Hence , 

x 2 “ Y 

As before, 


f(h. R \ x 2 

\ \ R 0 RoV 2 * 2 ) 


/X3X5 

\ *4 


1 “ ~ 


88 211 

/X4X5 


— )) 

*4*5 ) J 


AR 


x 4 * X 3 - — / X4X5 


(D16) 


(D17) 


Now 


U)* 


1 t * 
L '</ <PJ 


- P3) dr + 


(D18) 


p 3 - ^§7 (1 + 0.7746 <Ji PN) 


(D19) 


41 is different in reverse flow, 

♦ - KN $2 + 0.33 (D20) 

From (D13 ) , (D19), and (D20) 

p 3 " 0.837 V 2 * 2 ( l + 0,2556 PN2 + 0,7746 ™ " x^xf ' 2 ) (D21) 

From (D9), (Di8), and (D21) 


2 

• R X 2 X3 X5 

x 3 ■ L 0.837 V 2 ( M + “ 810 5 ' vj 

where 

89 - 1 + 0.2556 PN 2 


(D22) 


(D23) 


810 “ 0 . 77^6 g 2 

Now from 


(D24) 
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(WT) 3 - Y (d»2 T 3 - u 3 T 3 )dt + (WT) 3fl 
and Equation (B14), 



This concludes the derivation of the fifth order reverse flow 
system Equations. 

Summary of Reverse Flow Stalled Dynamics (RS) 



where 


D-4 


(D25) 

stalled 

(RSI) 

(RS2) 

(RS3) 

(RS4) 

(RS5) 
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corrected flow h r 


*3 
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APPENDIX E 




Figure E(a) Surge oscillations for PN»1 occur at A*26.96. 

(They also occur at A-27.18j they are not included here.) 


E-l 
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.at .it -a -a t>.» j * it .it 

nrc; (SfcoNBS) nrc isecokbj 


Figure E(c) Surge oscillations for PN».8 occur at A-30.309 
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Figure E(d) Surge oscillations for PN-.7 occur at A-29.072 
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